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Abstract: We have used diffraction gratings to simplify the fabrication, 

and dramatically increase the atomic collection efficiency, of magneto¬ 
optical traps using micro-fabricated optics. The atom number enhancement 
was mainly due to the increased beam capture volume, afforded by the 
large area (4cm 2 ) shallow etch 200nm) binary grating chips. Here 
we provide a detailed theoretical and experimental investigation of the 
on-chip magneto-optical trap temperature and density in four different 
chip geometries using 87 Rb, whilst studying effects due to MOT radiation 
pressure imbalance. With optimal initial MOTs on two of the chips we 
obtain both large atom number (2 x 10 7 ) and sub-Doppler temperatures 
(50 /iK) after optical molasses. 
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1. Introduction 

Laser cooling ||T]| revolutionised atomic physics by making it possible to rapidly and robustly 
chill dilute clouds of atoms to below a millionth of room temperature. At such slow speeds it is 
possible to observe and interact with atoms a thousand times longer. This makes measurements 
of atomic transitions correspondingly more accurate, suitable for advanced quantum metrolog¬ 
ical devices. The starting point for the majority of cold atom experiments is a magneto-optical 
trap (MOT) (2|[3j. A MOT both cools and traps atoms surrounded in red-detuned laser light 
via changes in the overall radiation pressure due to the Doppler and Zeeman effects, respec- 





tively. MOTs contain from one Q to « 10 11 [ 5,6 ] atoms at milliKelvin typical temperatures, 
and the atomic species cooled usually has a simple energy level structure - which covers all 
alkali metals, alkaline earth species and metastable noble gases. Other atomic species can also 
be very effectively cooled (e.g. holmium 0). and very recently the first molecular MOT was 
obtained ©• 

We have recently highlighted the benefits of using surface-patterned chips to make MOTs 
with relatively large numbers of ultracold atoms 0 , compared to prior micro-fabricated 
magneto-optical traps (MOTs) (K^[TT]]. The design requires only the alignment of a single 
input laser beam, and extends our previous work on tetrahedral mirror |T 2 ) and grating G3 
geometries, which are four-beam m equivalents of pyramidal MOTs p4) . Figure [I] depicts a 
‘like-for-like’ graphical comparison of the optical setup in a standard 6 -beam MOT (Fig.[lja) 
01313). a mirror MOT (Fig.Qb) 119]) and a grating MOT (GMOT, Fig.[I]c) ©). 

Here we study the atom-number characteristics of four different chip systems, and more com¬ 
plete details of different fabrication methods for the chips are discussed by Cotter et al. p0| . 
The two chip geometries we used are shown in Fig. [2] One chip style (Fig. |2ja)) comprises 
three one-dimensional gratings arranged in a pattern with 3-fold rotational symmetry (we used 
chips with periods of d = 1470nm and d = 1200nm). The other chip style (Fig. [2jb)) is a 
two-dimensional grating generated with circular regions centred in square units cells of side 
length d = 1080nm. Hereafter the chips will be referred to as TRI 15 , TRI 12 (Fig.[2ja)) and CIR 
(Fig.[2jb)). Although we have used chips TRI 15 and CIR previously (chips B and C in Nshii et 
al. (9j), purely for number characterisation, here we present density and temperature data for 
these chips. Moreover, we also investigate MOTs with improved capture efficiency for the ID 
structure (i.e. the new chip TRI 12 ) and improved radiation balance for the 2 D structure (chip 
CIRjqp): identical to CIR but with an input beam intensity spatially sculpted with a neutral 
density filter - for which the rationale is given in Sec. [4]). 

We begin by introducing simple theoretical models for the MOT temperature (Sec. [2]) and 
atom number (Sec. [3]), then consider the amelioration of any effects caused by imbalanced 



Fig. 1. Magneto-optical trap geometries: six-beam MOT (a), reflection MOT (b) and grat¬ 
ing MOT (c). The black cuboid indicates a vacuum cell, yellow fibers deliver cooling and 
repumping light (red) which is collimated by lenses (blue) and appropriately circularly- 
polarised with quarter waveplates (green). The vacuum pump, atom source, anti-Helmholtz 
magnetic coils etc. are omitted for clarity. 













Fig. 2. The ‘zoomed in’ microstructure on the different 2 x 2 cm 2 chips used in our studies. 
The black and white zones are separated by a height difference of ~ 200nm, one quarter 
of the wavelength of the rubidium cooling light (780nm). Chips TRJ 45 and TRI 12 have the 
pattern shown in (a), comprising three one-dimensional gratings with period d (blue line) 
of 1470 and 1200nm, respectively. The other chip, CIR, has the pattern shown in (b), with a 
grating unit cell (blue square) side length d = 1080nm. First order Bragg-diffracted beams 
have angles 32° (TRI 15 ), 41° (TRI 12 ), 46° (CIR) with respect to the chip surface normal. 


radiation pressure in grating MOTs (Sec. [4]). We then present the experimental data (Sec. [5]), 
comparing and contrasting with the theoretical predictions. Finally we present ‘best of both 
worlds’ findings (Sec. [6]) — by permitting analog temporal evolution of both MOT intensity 
and detuning we go beyond Nshii et al. |9| to achieve high atom numbers and sub-Doppler 
temperatures simultaneously using both ID and 2D grating structures. 


2. Theory: Doppler temperature 

Whilst sub-Doppler cold atom temperatures can be reached in optical molasses (Sec. HHH- 
[27) , the conditions for atom collection in a MOT are in stark contrast to those required to 
reach the lowest temperatures: a MOT needs high beam intensities relative to the saturation 
intensity (/ >> Is), detuning of a few natural linewidths A « — 2T and a magnetic gradient 
~ lOG/cm for a typical alkali metal vapour cell MOT with 25 mm beam diameter. For this 
reason the considerably simpler Doppler theory is a good description for our MOT temperature. 
The standard 3D Doppler result for a sample of atoms uniformly and isotopropically illuminated 
in three dimensions by light of total intensity Ij , with corresponding saturation parameter )3t = 
It/I s, and detuning A [21 ] is: 


Tu = T D 


l+jS r +4A 2 /r 2 

-4A/r 


(i) 


where Tp = In our 87 Rb experiment we cool using the rubidium D2 transition which has a 
natural linewidth T = 27T x 6.07 MHz. The standard Doppler temperature minimum in the limit 
of low intensity and detuning A = — T/2 is Jb, i.e. 146/iK for rubidium, with the more gen¬ 
eral lowest temperature 7b y / TT/3r minimised at a detuning A = — %/1 + fir T/2. Throughout 
this paper, for consistency, we use the saturation intensity averaged over all polarisations and 
magnetic sub-levels Is = 3.57mW/cm 2 - however we note that using the stretched state satura¬ 
tion intensity Is = 1.67 mW/cm 2 for the theoretical temperature of Eq. {!]) leads to much better 
agreement with our experimentally measured temperatures (Sec. [5]). 

Although the velocity distribution in a magneto-optical trap is Gaussian, the system is by 
no means in thermal equilibrium in the usual sense, with the typical light scattering rate from 
atoms vastly larger than the interatomic elastic collision rate. The final ‘temperature’ in a given 
direction is given by the balance between the heating (cf. diffusion coefficient and light scatte- 














ring rate) and the velocity damping constant y. In our system the effective temperature parallel 
to and perpendicular to the grating are therefore expected to be different. 

For our gratings if the incident beam has a spatially uniform intensity Iq and wavevector 
ko = k{0,0, —1}, for balanced optical pressure from N diffracted orders each should ideally 
carry an intensity of / z = Iq/(Ncos(X ) with i G {1, The wavevectors of the diffracted 

beams are k z = k{sin a sin , sin a cos , cos a}, where a is the angle the diffracted orders 
make to the grating normal {O,O,l}(a = 20cf. |l2j[l3]). The total intensity is therefore 7/ = 
Y 4 L 0 h = (1 + sec a )k ), with each beam of intensity 7 Z contributing a relative intensity-dependent 
heating proportional to its intensity along the beam direction due to absorption, and an equal 
associated heating rate due to spontaneous emission. The heating from spontaneous emission 
is assumed for simplicity to be isotropic, with the heating rate for light with a wavevector 
k z = {ki x , ki y , ki z } proportional to h 2 Rik 2 {^ in the x,y and z dimensions, for a scattering 
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Assuming dipolar (as opposed to isotropic) spontaneous emission from a circularly polarised 
atom adds complexity and does not significantly alter the predicted temperatures. The heating 
due to absorption is apportioned ti 2 Ri{ki 2 ,ki y 2 ,ki 2 } in the three spatial dimensions. The total 
heating rate near zero velocity in each dimension is thus: 
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which simplifies to 


D = h 2 k 2 ^-R ({2,2,8} + sec a {3 sin 2 a + 2,3 sin 2 a + 2,6cos 2 a - 
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In balanced laser cooling the relative damping constants can be determined by the Taylor 
expansion of the total force equation F = about v = 0. The damping constants in 

the xy directions (y| parallel to the grating) and z direction (y_j_ perpendicular to the grating) are 
thus given by m 
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{ 7 ||, 7 ±} 


—7-v = —{Y\\v x ,Y\\Vy,y_\_v z } 

f sin a tan a 1 + cos a 1 
761 -4-»- 2 -) ’ where % 


-16M 0 k 2 RA 

/ 5 r 2 (i + jS r + ^) 


( 5 ) 


is the damping force in a standard 6 -beam MOT. The temperature in each dimension can then 
be determined as proportional to the ratio of the diffusion D (Eq. |4])) to the damping y (Eq. ([5])): 

{7]|,71} = ^{3 + csc 2 (a/2),3 + seca}, ( 6 ) 


which reduces to the standard isotropic temperature Ty when a = arccos(l/3). Note this tem¬ 
perature only depends on the grating diffraction angle and not on the number of diffracted 
orders N. 

The temperatures obtained experimentally (Sec. [5]) do not reflect the theoretical temperature 
disparity in the directions parallel and perpendicular to the grating. The model above only 
includes forces due to absorption of light directly from the laser beams, however for MOTs with 
larger atom number and light scattering rate (like in our experiment), there will be a significant 
contribution to the forces in the MOT due to the l/r 2 repulsion from closely packed atoms with 






inter-particle spacing r due to absorption of spontaneously emitted light from neighbouring 
atoms | [28] - [3T| . These forces effectively mix the energies across dimensions, creating a more 
uniform temperature distribution. The theoretical temperature we use in both dimensions is 
therefore T = |7j| + ^7j_, which agrees surprisingly well with experimental values given the 
simplicity of the model. 

For completeness and to determine theoretical MOT density, one can also derive the spring 
constant K of the trap, expanding the total force F, but this time relative to position. The a 
dependencies of the damping constants are {fCy, k±} oc {i sin a tan a, 1 — cos a}, the same as 
derived by Vangeleyn et ah in [12]. The damping therefore reduces to an isotropic | (^) of 
the standard 6-beam MOT radial (axial) K values, if the GMOT has a = arccos(l/3) and laser 
beams of the same intensity as the standard MOT. Under these conditions the GMOT should 
be spatially isotropic - it will have equal rms cloud radii as a standard MOT in both directions 
parallel to the grating, but twice the standard MOT extent in the direction perpendicular to the 
grating. An important caution against using this Doppler spring constant to determine density, 
however, is that the experimental spring constant of a trap is notoriously difficult to determine, 
and the reradiative forces f28] - f3T| discussed in the previous paragraph make simple models of 
the effective spring constant considerably more difficult. 


3. Theory: Doppler atom number 

In order to predict the number of atoms in the MOT, we use the standard method of Lindquist 
et al. fl6| : treat the MOT as a spherical target region with radius r and cross-section 4;rr 2 , then 
find the flux of atoms incident on the target, with speeds less than the capture velocity v c of the 
MOT slowing region. In the ideal case of a rubidium-dominated vapour both loading and loss 
mechanisms are proportional to atomic density, which is therefore irrelevant for calculating 
atom number. It is necessary to know the collisional cross section of cold atoms as seen by 
hot atoms, which for rubidium is a = 2 x 10 -17 m 2 [32] (similar to caesium [16]). We make 
the assumption of a ID optical molasses from two counter-propagating laser beams. The total 
acceleration is modified by a prefactor rj to allow for geometric effects: 
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By rearranging one arrives at an equation of the form f(v)dv = Bdx , where f(v) is a quartic 
polynomial in v and B is a constant. This equation can then be integrated to analytically de¬ 
termine x(y). The integration constant is chosen such that x(0) = 0. Setting x(y c ) = 2 r (atoms 
are stopped over a distance of twice the MOT radius) and subsequently inverting gives v c , the 
capture velocity. The total steady-state number of atoms in the MOT is then fl6| : 


4k r 2 f v c 
8(7 \Vt 
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( 8 ) 


where vj = \j2ksT/m is a thermal velocity in the (hot) background vapour. 

We note that the saturation intensity features in the theoretical temperature and both theoret¬ 
ical and experimental atom number, and we therefore consistently use the saturation intensity 
averaged over ah polarisations and magnetic sub-levels Is = 3.57mW/cm 2 , which is valid for 
small detunings |23|24| . There are, however, arguments for using a value closer to the stretched 
state saturation intensity Is = 1.67mW/cm 2 [ 22J. We used this value in our previous work on 
atom number (9) as a precaution to prevent over-counting the atom number, although for high 
levels of saturation the effect of changing Is on experimental atom number is negligible. 











4. Theory: Non-uniform beam intensity 

In the theory so far we have assumed a spatially uniform incident laser intensity. However, 
real propagating laser beams tend to have Gaussian spatial intensity profiles. If such inten¬ 
sity distributions are apertured, or their intensity is spatially modified in other respects, as 
long as the modifications occur on large (centimetre size) beams, diffraction effects can be 
neglected. While various refractive beam shapers which transform Gaussian beams into flat- 
top intensity beams are available on the market, these tend to be expensive (albeit power ef¬ 
ficient) and designed for relatively small output beam diameters. Spatially shaping a Gaus¬ 
sian beam with radial profile I = /oexp(—2 r 2 /w 2 ) using an apodizing filter of outer radius 
r a with a neutral density ND = exp(2(r 2 — r a 2 )/w 2 ) leads to a uniform beam intensity, with 
reduced intensity /oexp(— 2r a 2 /w 2 ) across the beam. The relative power in the beam is then 
2(r a /w) 2 exp(—2 r a 2 /w 2 ) which has a maximum of \/e « 37% when the aperture radius is 
r a = w/y/l - i.e. even under optimal conditions it is a fairly lossy solution. 

To determine the effect of Gaussian intensity profiles, we fully integrate the acceleration from 
a MOT with a single input beam and 3 diffracted orders at diffraction angle 41°, just like Chip 
TRI 12 (Fig. [3}. We assume similar parameters to the MOT experiment - a Gaussian input trap 
beam with intensity 10mW/cm 2 and detuning A/T = —1.7 apertured to a 12 mm radius and a 
magnetic quadrupole with axial gradient 15 G/cm. In situations where the beam intensity is spa¬ 
tially uniform (Fig.[3ja)) the radiation pressure balance is uniform throughout the beam-overlap 
volume. However, with a Gaussian beam intensity profile (Fig. [3jb)), the diffracted orders at 
the MOT location originate from lower and lower intensity regions of the Gaussian beam as 
the MOT position is raised from the grating surface. This means that the MOT is pushed down¬ 
wards more, relative to the magnetic quadrupole centre, as the MOT is raised above the grating 
surface. Moreover, there is a marked change in the relative trapping and damping constants: 
{K||,Kl,y||,yx,} = {0.21,0.37,0.21,1.31} in Fig.^a) and {0.10,0.49,0.16,0.96} in Fig.pjb) 
(compared to the isotropic values of {1,1,1,1} for a = arccos(l/3) and uniform intensityf 

The case of uniform intensity is preferable, as long as one has optimal diffraction efficiency 
in the gratings, as then the MOT and molasses properties do not depend on MOT position. How- 
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Fig. 3. The full theoretical acceleration above chip TRI 12 , shown both as vector streamlines 
and magnitude. The white arrow indicates the centre and direction of the beam incident on 
the grating (gold line at figure base), whereas the white cross indicates the quadrupole 
magnetic field centre. In figure a, a uniform beam intensity profile is used, whereas in b we 
use a Gaussian beam with \/e 2 radius of 2 cm, like in the experiment. Black dashed lines 
indicate the boundaries of the beam overlap (capture) volume. 














ever, if one uses Gaussian beams and the diffraction efficiency is too high, one can compensate 
by raising the MOT position to a location where optical forces are balanced and suitable for 
sub-Doppler cooling. In Sections [5] and [ 6 ] with chip CIRjqp) we show how locally shaping the 
input beam intensity can dramatically improve both MOT and molasses performance. 

5. Experiment 

The experimental conditions were very similar to those detailed in the Methods section of 
Nshii et al. (9). One major difference is that we now have analogue as opposed to binary in¬ 
tensity control, using an acousto-optic modulator (AOM) after our tapered amplifier prior to 
fiber-coupling. After the MOT reached steady-state atom number, the atoms were ballistically 
released for 0, 1.5 or 2ms, and Gaussian fits were used to extrapolate cloud rms radii and 
thereby MOT temperatures. To build up statistics each drop time was iterated five times. For 
intensities above and below 10mW/cm 2 the fluorescent imaging exposure times were 430/is 
and 150/is, respectively. 

Figures [4] (Chips TRI 15 and TRI 12 ) and [5] (Chips CIR and CIR^p)) illustrate the intensity- 
and detuning-dependence of a variety of MOT parameters: theoretical atom number (a), ex¬ 
perimental MOT atomic number (b) and spatial density (c) as well as the theoretical Doppler 
temperature (d) with the temperatures experimentally measured in the directions both paral¬ 
lel (e) and perpendicular (f) to the grating. Note that the theoretical atom number saturates at 
about three times lower intensity than that observed experimentally - we seem to need three 
times more laser power experimentally than we would theoretically expect. This discrepancy 
may be considered relatively minor given the very simple nature of the theoretical model, how¬ 
ever we believe further investigation, and full comparison to a standard 6 -beam MOT may still 
be warranted. 

In prior measurements of T± on Chip CIR (Fig. [5jf) CIR) there was an even larger discrep¬ 
ancy with theory (very high experimental temperatures) due to imbalanced optical force during 
the fluorescence imaging. This was determined via the noticeable centre-of-mass velocity ac¬ 
crued in the direction perpendicular to the grating during time-of-flight. In all the new data 
for Figs. [4] and [5] any imbalanced radiation pressure (mainly for Fig. [ 5 ]) was balanced using a 
non-zero magnetic field oriented perpendicular to the grating. 

Note the significant increase in atom number (x3) and density (x 6 ), with lower atomic 
temperature on chip TRI 12 compared to chip TRI 15 - in future we intend to test gratings with 
even smaller period (larger diffraction angle) to determine if this trend continues. 

6. Experiment: best of both worlds 

For many metrological experiments, such as atomic clock measurements, the precision achieved 
strongly depends upon the temperature and total population of the atomic ensemble. Hence, an 
understanding of the final number of atoms that can be brought from the MOT stage to ultracold 
temperatures provides a good indication of the full capability of the GMOT. In Nshii et al. ||9) 
we effectively had only binary control of the intensity - this meant we could either get a large 
number of atoms on Chip CIR with high intensity, or good optical molasses on chip TRI 15 with 
lower intensity light (and a correspondingly small number of atoms loaded in the MOT). 

After analysis of the results presented in Section [5j the most suitable gratings for MOT 
number and temperature are Chip CIR^jy and TRI 12 . To study the number of atoms that 
could be brought to sub-Doppler temperatures, a MOT of intensity 40mW/cm 2 and detun¬ 
ing of A/r = —2.2 (initial MOT indicated by black dots in experimental Figs.[ 4 ]TRIi 2 and [ 5 ] 
CIR^jy) was transferred into optical molasses. The molasses consisted of turning off trapping 
magnetic fields, whilst using a two-step frequency jump of A/r = —5.3 and A/r = —8.3, each 
of duration 12 ms. 
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Fig. 4. The full GMOT characterisation when using chips TRI 15 and TRI 12 with beam 
diffraction angles a = 32° and 41°, respectively. The labels a-f indicate theoretical atom 
number (a), experimental atom number (b) and spatial density (c), the theoretical tempera¬ 
ture (d) and experimental temperatures in directions parallel (e) and perpendicular (f) to the 
grating. Light gray indicates indeterminate or out-of-range data. Medial shows the theo¬ 
retical temperature (d) for 1$ = 1.67mW/cm 2 , which has clearer experimental agreement. 
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Fig. 5. The full GMOT detuning and intensity characterisation when using chip CIR (stan¬ 
dard Gaussian input beam a-f) and CIRj^p) (CIR with a T = 0.74 ND filter in the cen¬ 
tral 3 mm diameter of the input beam, al-fl), both with beam diffraction angle a = 46°. 
Note the slight increase in atom number, and significant decrease in temperature using chip 
CIRnd relative to chip CIR. Media2 is the Is = 1.67mW/cm 2 version of the theoretical 
temperature (images (d)). 





































































Fig. 6. Temperature measurements in directions both parallel (red) and perpendicular (blue) 
to the grating for chip TRI 12 (a) and CIR^d (b). Both chips give a 3D average temperature 
T = |7|| + ir ± = 50MK. 


Under these conditions, Chip CIR^pj brought 2 x 10 7 atoms down to a 3D average tem¬ 
perature of 50 fl K (Fig. | 6 ja)). Chip TRI 12 used the same settings with the exception that the 
molasses intensity was reduced to 23mW/cm 2 - bringing 2.5 x 10 7 atoms down to a 3D average 
temperature of 50 /iK (Fig.^b)). Both temperatures are similar to those achieved by Nshii et 
al. in j9| (with TRI 15 ), but with an order of magnitude more atoms after molasses (N ~ 2 x 10 6 
in [9), due to the lack of intensity control discussed above). Temperature T is determined by 
fitting a 2 = Gq 2 +ksTt 2 /m to the data, where G is the standard deviation of the Gaussian fit to 
the atomic cloud spatial distribution, m is the mass of an 87 Rb atom, ks is Boltzmann’s constant 
and do is the initial (t = 0 ) standard deviation of the cloud. 

7. Interpretation and conclusions 

Both the experimental MOT atom number and cloud temperature demonstrate surprisingly 
good agreement with theory, given the simplicity of the theoretical models employed. 

We note that higher densities are achieved on Chip CIR and CIR^jy, largely due to the 
larger diffraction angle a, with higher associated damping 7 and trapping K coefficients. There 
is a dramatic density improvement with TRI^ over TRI^ and this trend may continue with 
higher values of a toward the tetrahedral configuration a = arccos(l/3), with the caveat that 
the capture volume and grating diffraction efficiency will decrease. In Doppler theory the cloud 
size scales as g oc \/YT /k, however we note from the relatively constant experimental MOT 
density for a variety of detunings and intensities that we are reaching sufficient atom numbers 
to likely be entering the MOT ‘constant density regime’ where reradiation forces become more 
important than restoring forces |28j-[3i~[ . 

In terms of the temperature, significant sub-Doppler effects are not present, however we note 
that one would only expect sub-Doppler cooling under conditions of low intensity and large 
negative detunings, where we have insufficient data coverage to draw any further conclusions. 
The atom number is sufficiently large that effects where the MOT temperature is expected to 
increase with atom number (T oc N 1 / 2 ) might also become evident [33 34]. 

Compact quantum measurement devices are a burgeoning area (35 J37|, with grating tech¬ 
nologies providing a new way to optically simplify transportable ultracold atom experiments. 
In future we hope to investigate the effect of grating diffraction angle a on both grating optical 
properties (38) and cooling performance. The current work could also be extended from inves¬ 
tigating MOTs into sub-Doppler optical molasses regimes (39) , preparing the path for future 
studies of on-chip Bose-Einstein condensates and extremely stable optical lattices (9). 
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